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The design of an electrodeless thruster that relies on a pulsed, rf-assisted discharge 
and electromagnetic acceleration using an inductive coil is presented. The thruster design 
is optimized using known performance sealing parameters and experimentally-determined 
design rules, with design targets for discharge energy, plasma exhaust velocity, and thrust 
efficiency of 100 J/pulse, 25 km/s, and 50%, respectively. Propellant is injected using a 
high-speed gas valve and preionized by a pulsed- RF signal supplied by a vector inversion 
generator, allowing for current sheet formation at lower discharge voltages and energies 
relative to pulsed inductive accelerators that do not employ preionization. The acceleration 
coil is designed to possess an inductance of at least 700 nH while the target stray (non- 
coil) inductance in the circuit is 70 nH. A Bernardes and Merryman pulsed power train 
or a pulse compression power train provide current to the acceleration coil and solid-state 
components are used to switch both powertrains. 




Nomenclature 

A c 

magnetic core area, m 1 2 

m bi t 

mass bit, kg/pulse 

B s 

saturation magnetic field, T 

n 

number of turns 

C 

capacitance, F 

OD 

outer diameter, m 

I 

current, A 

t ' 

time, s 

ID 

inner diameter, m 

V 

voltage, V 

3 

linear current density, A/m 

Vo 

initial charge voltage, V 

Z . 

coil width, m 

v z 

exhaust velocity, m/s 

L 

dynamic impedance, H/s 

z 

axial distance, m 

L' 

inductance per unit length, H/m 

ZQ 

electromagnetic decoupling distance, m 

L* 

inductance ratio 

a 

dynamic impedance parameter 

Lo 

L C 

stray inductance, H 
coil inductance, H 

Vt 

thrust efficiency 
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L Introduction 


P ulsed inductive plasma accelerators are spacecraft propulsion devices in which energy is stored in a 
capacitor and then discharged through an inductive coil. The device is electrodeless, inducing a current 
in a plasma located near the face of the coil. The propellant is accelerated and expelled at a high exhaust 
velocity (O (10 km/s)) by the Lorentz body force arising from the interaction of the plasma current and the 
induced magnetic field. Presently, there are two concepts which operate on this principle. One is the Pulsed 
Inductive Thruster (PIT), 1 in which both propellant ionization and acceleration are performed by the pulse 
of current flowing through the inductive coil. In contrast, the Faraday Accelerator with Radio-frequency 
Assisted Discharge (FARAD) 2 * uses a separate inductive discharge to preionize the propellant before it is 
accelerated by the current pulse in the coil. 

Inductive plasma accelerators are attractive as propulsive devices for many reasons. The lifetime and 
contamination issues associated with electrode erosion in conventional pulsed plasma thrusters (PPTs) do 
not exist in devices where the discharge is inductively driven. In addition, a wider variety of propellants 
(e.g. CO 2 , H 2 O) becomes available for- use when compatibility with metallic electrodes is no longer an 
issue. Moreover, pulsed inductive accelerators (indeed, pulsed accelerators in general) can maintain the 
same performance level over a wide range of input power levels by adjusting the pulse rate. 

The state-of-the-art PIT thruster is 1 m in 'diameter and requires approximately 4 k J/pulse for efficient 
operation (using ammonia as a propellant). 1 In addition, this thruster has multiple capacitors which must 
be individually (and simultaneously) discharged to produce the azimuthal current in the coil. Recently, a 
FARAD-based benchtop accelerator was used to demonstrate that preionization of the propellant using a 
steady-state helicon discharge allows for formation and. acceleration of an inductive current sheet at much 
lower discharge energies and voltages than those found in the PIT. 2,3 While the experimental apparatus in 
that' study was not optimized as a thruster, the work resulted in the analytical derivation of performance 
scaling parameters and the formulation of a set of design rules for a FARAD thruster. 3-5 % 

In this paper, we present the design of a complete laboratory- model FARAD thruster with all the sub- 
systems (mass injection, preionization, and acceleration) integrated into a single unit. The system uses a 
pulsed, ringing RF- frequency signal produced by a discharging vector inversion generator (VIG) to preionize 
the gas. The acceleration stage operates on the order of 100 J/pulse and can be driven by several different 
pulsed powertrains. These include a simple capacitor coupled to the system, a Bernardes and Merryman 
based drive circuit, and a pulse-compression circuit that takes a temporally broad, low-current pulse and 
transforms it into a short, high-current pulse. 

The outline for the rest of this paper is as follows. In Sect. II we review the scaling parameters and 
design rules for a FARAD thruster and show how we arrived at our design. The designs of the individual 
subsystems in our FARAD thruster are discussed in Sect. III. Any data obtained or modeling performed to 
test certain subsystems and validate our design approach are included in this section as well. 

II. Scaling and Design Considerations 

The inductive coil and its associated drive circuit are optimized for plasma acceleration and current 
sheet formation. From the very beginning, we assumed a discharge energy of 100 J/pulse and fixed our 
acceleration coil dimensions with an inner diameter of 12 cm and an outer diameter of 30 cm. We proceed 
with a discussion of the performance scaling parameters and the optimization procedure. 


A. Acceleration Optimization 


Optimization of the acceleration process follows directly from a set of recently derived performance scaling 
parameters. 4 Specifically, we focus on the dynamic impedance parameter, a , which is defined in terms of 
measurable quantities as 


C?V§Lc 
2m b it-?o ’ 



where C is the capacitance in the bank, Vo is the charge voltage, Lc is the inductance of the coil, mbit is 
the mass bit or propellant mass expended in a single pulse, and z 0 is the decoupling distance or length scale 
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over which inductive plasma acceleration occurs. The parameter a can be recast as 


1 CV?/2 L , f 2nVLoC \ 2 

8?r 2 * m hi tV%/2 V L 0 /L ) ’ 


(2) 


where v z is the propellant exhaust velocity, Lq is the stray (non-coil) inductance in the circuit, L* is the 
inductance ratio Lq/Lc, and X is the dynamic impedance which is defined as v z U where L' is an effective 
inductance per unit length equal to Lq/zq. We recognize several important ratios in Eq. (2), specifically the 
inverse of thrust efficiency (defined as the ratio of directed jet kinetic energy to input electrical energy); the 
Lovberg ratio, 6 which is !;/£*; and the ratio of the resonant period of the unloaded circuit (i.e. without the 
coil inductance) to the time scale that the current sheet remains in the acceleration region before decoupling 
from the coil. The last ratio can also be understood as a quantitative measure of the transfer of stored 
electrical energy into directed kinetic energy. ' 

Previous work has shown 4 that in an efficient pulsed inductive accelerator a should be between 1 and 3. 
We choose 3 as a target value for our calculations based on the expected mass distribution in front of the coil. 
In addition, we set a target exhaust velocity of 25 km/s and a thrust efficiency r} t of 50%. For an efficient 
pulsed electromagnetic accelerator, the value of X* should be less than unity. We felt that a value of 0.1 was 
realizable in the current configuration, and based our analysis on Lq — 700 nil and Lq = 70 nH. Finally, 
measurements indicate that for our coil dimensions, the electromagnetic coupling distance zo is roughly 3.7 
cm. Using these numbers, we find that the unloaded circuit natural period is 5.1 ps and the capacitance 
required to given this period is 9.4 jjlF. Our capacitor bank is a modular design and each capacitor has a 
value of 2 /xF, so we choose our target value C = 10 /jF. 


B. Current Sheet Formation 

When the current pulse in the induction coil is initiated, it must transform the preionized plasma into a 
highly ionized, magnetically impermeable current sheet with essentially no delay. If this does not happen, 
.the magnetic field produced by the induction coil will rapidly diffuse through the highly resistive plasma and 
radiate away, performing little or no work on the propellant. It was recently proposed 3, 5 that formation of 
a current sheet capable of supporting efficient acceleration will occur if the linear current density rise rate 
in the coil is roughly 10 12 A/(m s). If the coil is electrically connected to a single capacitor bank and is 
comprised of multiple spirals each making n turns, then the initial current rise rate required from the bank 
is 

dJ _Cdj_ 
dt t=0 . n dt 

where dj/dt is the linear current density rise rate and £ is equal to the distance from the inner to the outer 
radius of the coil (the length over which the linear current density is defined). Following this argument, the 
powertrain must produce an initial current rise rate equal to 45 kA//xs. If we assume the acceleration coil 
is being driven by a powertrain approximated by an RLC circuit, then the initial current rise rate is Vo/L 0 . 
Furthermore, if T 0 — 70 nH, the initial capacitor charge voltage Vo should be at least 3150 V. 

We summarize our target design parameters in Table 1 . 



Table 1. Summary of target FARAD thruster design parameters. 


V z 

25 km/s 

Lq 

70 nH 

Vt 

50% 

L c 

700 nH 

Energy 

100 J/pulse 

Zo 

3.7 cm 

a 

3 

c 

10 fjF 

Coil ID 
Coil OD 

12 cm 
30 cm 

Eo 

>3150 V 
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III. Thruster Design 


In the present FARAD thruster iteration, there are three major subsystems: propellant injection, preion- 
zation, and acceleration. Each of these are discussed in turn below. Unless otherwise noted, these components 
are being fabricated or integrated by Radiance Technologies. 

A. Propellant Injection 

The propellant injection system consists of two main components. These are a high speed gas valve from 
Applied Pulsed Power to pulse the gas into the system and a manifold designed to direct the radial flow 
outward across the surface of the thruster and distribute it as evenly as possible. The estimated distribution 
time for argon gas across the drive coil is 0.24 ms. 

The gas injection system uses a high speed valve (see Fig. 1A) with a radial orifice capable of generating 
gas pulses as short as 70 ^s in pulse width. This valve employs a low mass beryllium copper Belleville spring 
diaphragm that lifts off the seal at its outer diameter permitting gas to flow radially from the valve when 
open. The diaphragm is driven by a pancake solenoid, which is powered by solid-state switched 10 kA, 3 kV 
pulse generator. Continuous mode operation at 5 pulses per second or burst-pulse mode operation at 100 
Hz can be achieved without forced cooling. The expected lifetime of the valve is 10 8 pulses. 



Figure 1. A) High speed valve. B) Rendering of the gas distribution manifold. 


The manifold, seen in Fig. IB, is fabricated from Maeor and consists of 72 semicircular channels each 
0.127 mm in radius. The manifold material was chosen to support machining of the small gas distribution 
channels. Macor has the added benefit of having a high temperature range and hardness, both of which 
should make it more resistant to degradation in the presence of the ionized propellant. The combined 
cross sectional area of the channels was chosen to set the mass flow rate to roughly 625 mg/s based upon 
distributing 150 fig of argon gas over the drive coil in 0.24 ms. 

When the valye is opened, gas expands radially into the manifold channels. The channels are short 
enough relative to their diameter to ensure that wall friction has a minimal effect on the mass flow rate. The 
gas in the manifold does not begin to supersonically expand until it reaches the manifold exit where the Mach 
number is unity. Restricting the gas flow until it is at a larger radius (namely the ID of the acceleration coil) 
keeps the propellant closer to the drive coil over the entire radial range. Limiting the flow in the manifold 
to subsonic speeds also increases the distribution time in comparison to a completely supersonic expansion, 
resulting in more time between the initiation of the gas pulse and the other pulses in the system. 

Once the gas distribution system has been'built, a fast ionization gauge will be used to measure the time- 
varying propellant distribution over the acceleration coil face. This measurement will also help in estimating 
and adjusting the mass injected during each pulse. 
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B. Preionization 

An antenna assembly is used to capacitively couple energy into the propellant for the purpose of creating the 
preionized plasma. The antenna is connected to a vector inversion generator (VIG), which supplies energy in 
the form of a short, pulsed-RF signal. 7,8 Gas ionization experiments have shown that a significant portion 
of the propellant can he ionized for up to 20 /is using a VIG operating at as little as 30 mJ of input energy. 

The elements of the VIG are shown schematically in Fig. 2 A and B, and a VIG potted in castrol oil 
is shown in Fig. 2C. The lines shown in the schematics are metallic foils and the space between them is 
dielectric film. In this spiral configuration, the VIG consists of two parallel plate transmission lines wound 
on a mandrel and sharing a common conductor. During operation, the array is dc charged to approximately 
2 kV with energy stored in the capacitance of the spiral windings, C v i g . The erection sequence is initiated 
by closing the input switch, Si- ' 



Figure 2. Schematic of a spiral line VIG A) DC charged and B) fully erected. Solid lines are metal foils and 
the intermediate spaces are insulators. Arrows represent the electric field vectors in each configuration. C) A 
VIG potted in castrol insulating oil. 


When switch Si closes, an electromagnetic wave begins traveling up one of the strip lines. As a first- 
order approximation, due to the inductance disparity between the two conductors, a wave does not propagate 
quickly in the second line. The net result is that the VIG is comprised of two LC circuits - one “fast” and 
one “slow” ... The electric field vectors plotted in Fig. 2a transition to the fully-erected configuration shown in 
Fig. 2b after the “fast” wave travels to the “open end” and back to switch Si. If the transit time in the fast 
circuit is sufficiently short when compared to the slow circuit, the voltage across the array will be close to 
the ideal output of 2 nVo> where Vq is the initial charge voltage and n is the number of turns in the VIG. The 
net result is a transient high voltage between the lines at Si and the outermost turn on the VIG. The output 
voltage waveform is a ramp function with rise-time equal to the transit time up and back in the active line. 
Voltage outputs as great as 95% of ideal have been realized in certain geometries. 

C. Acceleration 

1. Drive Coil Design 

The coil design was performed at NASA-MSFC and is a critical element in fixing the overall performance 
of the FARAD thruster. A geometry consisting of several interwoven spiral elements is chosen as this yields 
a current density over the coil face that is approximately uniform. This leads to a relatively uniform radial 
magnetic field between the acceleration coil and the induced current sheet. The radially uniform field results 
in approximately equal axial acceleration profiles at any radial location in the current sheet, leading to full 
utilization of the propellant mass injected over the coil. To approach the performance levels assumed in 
Sect. II, the coil’s inductance and acceleration stroke length must be close to those given in Table 1. The 
coil inductance is important when one considers that the ratio of the coil inductance to the stray (non-coil) 
inductance in the circuit is representative of the amount of energy that can be electromagnetically deposited 
into acceleration of the plasma. 
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We surveyed several inductance approximation techniques in an 
attempt to find one that was appropriate for evaluating the induc- 
tance of a spiral acceleration coil like the one found in FARAD. These 
techniques included the Wheeler and modified Wheeler inductance 
formulas, 9 an approximation based on an azimuthal current sheet 
model, the use of conformal modules, 10 an approximation based on 
a flat spiral coil, 11 and an expression based on a curve fit of exper- 
imental data. The approximation techniques were able to predict 
the inductance of a single coil with relative precision. Unfortunately, 
they were not capable of modeling the inductance of a multiple- wire 
spiral coil possessing the dimensions of the FARAD thruster because 
of the addition of mutual inductance to the problem. 

Because of low confidence in the inductance approximation tech- 
niques for the FARAD coil, we performed a series of experimental 
inductance measurements using the coil shown in Fig. 3. The coil 
possessed multiple two turn spirals with the first turn of each coil 
traversing from the outer to inner diameter and the second turn 
returning to the beginning of the first turn. The two turns were spa- 
tially located on different planes to permit the wires to cross without 
difficulty. In this setup, we could test any number of coils from a 
single two-turn spiral to every spiral in parallel. 

We found that for two-turn spirals possessing an OD of 30 cm and 
an ID of 12 cm, the agreement between the experimentally measured 
coil inductance (that is, the inductance measured when the plasma 



Figure 3. FARAD test coil comprised 
of 12 two-turn coils and having an OD 
of 30 CM and an ID of 12 cm. 


is at £ =, oo) and that predicted using a non-coupled theoretical treatment (see Fig. 4A) is generally poor. 


This comparison is representative of the poor agreement between our measurements and all the theoretical 


models. 


From Figs. 4 A and B, we see that the measured inductance is roughly constant for 3 or more spirals. 
Physically, this result implies that 3 two- turn spirals in our configuration produces a magnetic field topology 
that reasonably approximates one produced by a purely azimuthal, uniform current distribution over the coil 
face. We elected to build a coil comprised of 6 two-turn spirals for the final design to remove nonuniformities 
in the near-field of the coil. The measured coil inductance in this configuration was 719 nH. 



0 1 2 3 4 5 6 7 0 2 4 6 8 10 12 

Number of Parallel Coils Number of Parallel Coils 


Figure 4. A) Comparison between the measured inductance (normalized) of a set of parallel, spiral coils and 
the predictions of a theoretical treatment assuming that the coils are uncoupled. B) Measured inductance 
(normalized) of a spiral coil as a function of the number of parallel coils used. Dashed lines inserted to 
demonstrate constancy of L co ji for a greater number of coils. 

. In the actual thruster, the inductance the coil presents to the rest of the circuit varies as the plasma 
moves away from the coil. We measured the coil inductance as a function of plasma position by placing 
a conducting plate at different axial positions relative to the coil to approximate the effect of the plasma. 
Inductance was measured using an Agilent 4285 A precision LCR meter, which is capable of measuring 
inductance at frequencies from 75 kHz to 30 MHz. Our data was collected primarily at 100 kHz as this is 
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close to the expected frequency of the current pulse in the FARAD. We found that the inductance presented 
to the circuit as a function of plasma axial position followed the well-known empirical scaling 12 

L (z) = L 0 + L c [1 - exp (-z/z Q )] , (4) 

where L is the total inductance measured in the circuit, Lq is the stray inductance in the circuit, and zq is 
the electromagnetic decoupling length scale. In our final design, this length scale was 3.7 cm, and in general 
our studies found that zq was always between 10 and 15% of the coil OD. A rendering of the final coil design 
is shown in Fig. 5. 

In addition to our inductance measurements, we 
modeled the coil using the 2-D axisymmetric mag- 
netostatic solver package in QuickField. We mod- 
eled the coil by taking a cross-sectional cut across 
the coil face and placing a circular conductor in ev- 
ery position in the r-z plane where our cut inter- 
sected the coil. This is not strictly correct since the 
coil is not axisymmetric. However, our approxima- 
tion should be valid so long as the magnetic field 
produced by the coil is close to that produced by 
a purely azimuthal, uniform current distribution at 
the coil face. 

We find that for a two-turn, 6 coil geometry with 
the cut taken at a location where the coil windings 
cross each other, the computed inductance from our 
model is 882 nH. It is not. unexpected that this num- 
ber is greater than the one measured experimentally. 

The coil can be completely isolated in the model 
and its inductance calculated, but. in the measure- 
ment the coil inductance is taken as the difference Figure 5. Rendering of the acceleration coil (shown 
between the inductance when the metal plate is at un P otted )* 

^ = .oo and when it is as close to the coil as physi- 
cally possible. A small amount of the coil inductance, representing of the difference between the metal plate 
being as close as possible and the plate actually being colocated with the coil, is inaccessible for performing 
electromagnetic acceleration on the plasma. In the experimental results, this inaccessible coil inductance 
is accounted for as additional stray inductance. There are also physical irregularities in the experimental 
setup where the coils cross each other and where they connect to coaxial feeds at the OD of the coil. A 
more refined assembly technique in the final FARAD acceleration coil should help to minimize the stray 
inductance attributable to these irregularities. 

2. Bernardes and Merryman Power Train 

The Bernades and- Merryman (B&M) circuit 13 is a pulsed power train designed to allow for energy recapture. 
The basic B&M circuit shown schematically in Fig. 6 A has two capacitor banks of equal capacitance. When 
the switch is closed, charge will oscillate between the two capacitor banks through the drive coil, represented 
in the schematic as Nl. At the end of the transient the voltage and energy on capacitor banks Ci and C 2 
will be equal. In this scheme, capacitor bank Ci can be brought back to full charge for the next pulse. 
Meanwhile, bank C 2 must be discharged before the next pulse is initiated. The excess energy on bank C 2 
could possibly to be used to power other spacecraft loads. This setup allows for recovery of at least half 
the energy remaining after the pulse transient. The target capacitance of 10 jjF required for acceleration 
optimization (see Sect. II) is given by the series configuration of the 20 /dF comprising both banks Ci and 
C 2 . However, the energy stored for each pulse is equal to l/2(Ci)k^ 2 . A target energy of 100 J/pulse is 
obtained for a charge voltage of 3160 V, which is slightly greater than the voltage required for current sheet 
formation. 

Example current and voltage waveforms for the B&M circuit are presented in Fig. 6B. These measure- 
ments demonstrate an additional advantage of the B&M circuit, namely that it precludes voltage reversal 
on either capacitor bank. This reduces the stress on the capacitors and increases their operational lifetimes. 
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Figure 6. A) B&M circuit schematic for the FARAD thruster. B) B&M voltage and current waveforms into 
a 540 nH inductive load for 1, 2, and 3 kV initial capacitor charge. 


To improve thruster lifetime relative to other pulsed inductive '■■ -****^ '*■* ' 

concepts that have been switched using spark gaps, 1 a solid-state, 
switch manufactured by Applied Pulsed Power was selected for in- 
tegration into the FARAD B&M circuit. To date, it has proved 
capable of switching 3 kV applied to bank Ci in the B&M circuit 
shown in Fig. 6 and discharged through a 540 nH inductive load 
with a peak current of 13 kA (see Fig. 6B). In the final FARAD 
thruster, the switch must be capable of switching a 3.2 kV capacitor j 
bank charge with a peak current over 30 kA. It must also possess low ; 
inductance so as to not place the thruster above the design target of j 
70 nH stray inductance. The capabilities of the switch have not yet 
been evaluated with respect to the FARAD thruster requirements. ! 

To further improve upon the B&M circuit’s energy recovery po- 

, , . i i . i . i . i i ■ , j. i a ■ u x« Figure 7. Applied Pulsed Power solid- 

tential, a latching design is also being studied. A circuit schematic st * te sw jtch 

of this concept in shown in Fig. 8A, where latching diodes have been 

added in parallel to two solid-state switches; one for each capacitor 

bank. It is well-known that plasma acceleration in a pulsed inductive plasma accelerator is most effectively 
accomplished during the first half-cycle, where coil current and dl/dt levels are greatest. Beyond that, a 
significant amount of energy is wasted in resistive heating of the external circuit and production of a mag- 
netic field that radiates into space. The diodes in Fig. 8 A serve to reduce the amount of energy wasted in 
the discharge by allowing current to flow in the circuit for only the first half-cycle of the discharge. During 
the discharge, the initially uncharged capacitor bank is partially charged by the current flowing through the 
coil from Gi'to C 2 . This reduces the amount of charge we must apply to bank C 2 before the next pulse, 
effectively recapturing energy that would otherwise have been lost in the system. 

A diode-latched design was bench tested with an inductive load equivalent to approximately 570 nH. 
The voltage and current waveforms from this test are shown in Fig. 8B. For a non-latching B&M circuit, 
the remaining energy would be equally split between the capacitor banks. In. the latching B&M test circuit, 
a larger percentage of the remaining energy in the circuit was recaptured by the second, initially uncharged 
capacitor bank. Testing is currently in progress to determine if this alternative design can be used at the 
charge, voltage and peak current levels required in the FARAD thruster. 


3. Pulse Compression Power Train 

A pulse compression system developed by the West Virginia High Technology Consortium Foundation pro- 
vides an attractive alternative means of powering the FARAD thruster. In this power train, initially slow, 
low-current pulses are compressed to yield short, high-current pulses. Solid-state switching is employed 
prior to pulse compression to keep current levels relatively low at the primary switching element. Typical 
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Figure 8, A) Diode-latched B&M circuit schematic for the FARAD thruster. B) Diode-latched B&:M voltage 
and current waveforms for 1, 2, and 3 kV initial capacitor charge. 


pre-compressed currents are low enough to allow the use of fast IGBT switches. The high pulse-repetition 
rates possible with IGBT switches lead to attractive values for thruster specific mass (kg/kw). 

A PSPICE model for the pulse compression circuit considered for this project is shown in Fig. 9A. Current 
from the primary energy-storage capacitor Cl is compressed a factor of five prior to application at the output 
capacitor C3. The resulting fast charge of capacitor C3 leads to low volt-seconds at the output magnetic 
switch L3. Low volt-seconds keep the core size for the magnetic switch reasonably small, as demonstrated 
using the integral form of Faraday’s law, 



Vdt, 


( 5 ) 


where A c is the area of the magnetic switch core, n is the number of turns made by the conductor wrapped 
around the core, V is the voltage applied to the conductor, and B s is the saturation magnetic field in the 
core. The magnetic switch can handle current amplitudes and rates-of-change of current well above levels 
tolerated by traditional solid-state switches. The drive coil for the thruster is represented in the model by 
the inductor L Coil and the resistor R3. The circuit inductance can be approximately transformed from the 
spatially-dependent form given in Eq. (4) to a temporally- varying form by using the substitutions z = v z t 
and zq ~ v z t 0 . This yields 

L(t) = Lo + Lo[l-exp(-t/t 0 )], - (6) 

where Lq is the stray drive coil inductance at t — 0, Lo + Lc is the total drive coil inductance in vacuum, and 
to is the characteristic time for plasma acceleration in the field of the drive coil. The time-dependent drive- 
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coil inductance is included in the simulation using the analog behavior models (not shown in the schematic) 
within PSP ICE. 

The circuit depicted in Fig. 9 A uses a novel energy-recovery scheme in which current passing through the 
drive coil is circulated back into the initial energy storage capacitor Cl through a resonant transfer from C3 
to Cl. Energy is held in Cl by means of steering diode D2. In effect, the circuit forms a pulse compression 
ring topology in which energy flows in only one direction around the circuit. This topology eliminates voltage 
reversal on the capacitive elements, resulting in a substantial increase in capacitor lifetime. Furthermore, 
recovered energy does not have to make a second transit of the circuit, as in other pulse compressors. 14 This 
feature leads to improved energy efficiency. 

Guided by PSPICE simulations, we are presently assembling the pulse compression ring shown schemat- 
ically in Fig. 9B. The, solid-state switching element in the circuit is’ formed from a circular array of eight 
IGBTs. The magnetic assist L10 is composed of eight independent windings on a single magnetic core. 
Each winding is connected to an individual IGBT. This construction assures uniform turn-on of all eight 
IGBTs and uniform current sharing between IGBTs. Steering diode D2 is formed from an array of six 
fast-response power diodes. The output magnetic switch L3 acts as a magnetic assist for this diode array to 
assure simultaneous turn-on of the array, low turn-on losses, and uniform current-sharing between individual 
diodes. 

IV. Design Summary 

We have presented the design of a complete FARAD thruster that is optimized based on previously derived 
performance scaling parameters and experimentally determined design rules. The thruster is designed to 
operate at a discharge energy of 100 J / pulse, and we have performed our design analysis using a target exhaust 
velocity and thrust efficiency of 25 km/s and 50%, respectively. Propellant is supplied to the thruster using 
a multi-channel gas distribution manifold that is fed by a high-speed, pulsed gas valve. Preionization is 
accomplished using a pulsed-RF signal supplied by a VIG. The coil has an OD of 30 cm and an ID of 12 
cm and is comprised of 4 two- turn spirals, with each conductor spiraling in on the first turn and back out 
on the second turn. Measurements and finite-element modeling have shown that a coil of this design will 
meet the target inductance of at least 700 nH. For the coil design, we have set our target stray (non-coil) 
inductance in the circuit equal to 70 nH to support both efficient acceleration and transformation of the 
preionized plasma into a highly-ionized, magnetically impermeable current sheet at t 0. The current is 
driven through the acceleration coil using either a Bernardes and Merryman pulsed power train or a pulse 
compression power train. Both power trains are switched using solid-state components, precluding the need 
for spark gap switches that are prone to erosion over time. 
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